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The spectrum and decay properties of radially excited D s states are examined in a new model. 
Good agreement is obtained with the properties of two recently announced D 3 mesons identified as 
■D s0 (2860) = cs(2P) and D*(2690) = cs as a possible mixture of (25; 3 Si) and (ID; 3 L>i). Searching 
for these mesons in B decays is advocated due to large predicted branching ratios. 

I. INTRODUCTION 

BaBar have recently announced the discovery of a new D s state seen in e + e~ collisions decaying to K~ir + K + , 
K-tt + tt°K+ (D°K+), or D+K^ l}. The Breit-Wigner mass of the new state is 

M(Z> sJ (2860)) = 2856.6 ± 1.5 ± 5.0 Mev (1) 

and the width is 

r(£> sJ (2860)) = 48 ± 7 ± 10 MeV. (2) 

The signal has a significance greater than 5 a in the D° channels and 2.8 a in the D + channel. There is no evidence 
of the £) sJ (2860) in the D*K decay mode 1 or the D s r] mode0. 

There is, furthermore, structure in the DK channel near 2700 MeV that yields Breit-Wigner parameters of 

M(L> sJ (2690)) = 2688 ± 4 ± 2 MeV (3) 

and 

r(D s j(2690)) = 112 ± 7 ± 36 MeV. (4) 

The significance of the signal was not stated. 

The discovery of these states is particularly germane to the structure of the D s (2317). For example, the low mass 
and isospin violating decay mode, D s ir°, of the I? s (2317) imply that the state could be a DK molecule 3J. If this is 
the case, the D s j(2690) could be a supernumerary scalar cs state. Alternatively, the D s (2317) could be the ground 
state scalar cs state and the new D s j's could be canonical radial excitations. Clearly, constructing a viable global 
model of all the D s states is important to developing a solid understanding of this enigmatic sector 

Previous efforts to understand the new BaBar states have argued that the D s j(2860) is a scalar cs state predicted 
at 2850 MeV in a coupled channel model 5] or that it is a J p = 3 _ cs stateQ. 

Here we pursue a simple model that assumes that all of the known D s states are dominated by simple cs quark 
content. It is known that this is difficult to achieve in the 'standard' constituent quark model with 0(a s ) spin- 
dependent mass shifts because the _D s o(2317) is much lighter than typical predictions (for example, Godfrey and Isgur 
obtain a D s o mass of 2480 MeV 7] ) . An essential feature in such phenomenology has been the assumption of two 
static potentials: a Lorentz scalar confining potential and a short range Coulombic vector potential. Following the 
discovery of the D s (2317), Cahn and Jackson ;8| analysed the D s states with a scalar potential S, whose shape they 
allowed to be arbitrary, while retaining a vector potential V that they assumed to be Coulombic. In the limit that 
the mass mi 3> m\ this enabled the spin dependent potential applicable to P-states to take the form 



V SD =XL-Si + 4tL-S 2 + tS u (5) 

(see the discussion around Eq. 1 of Q for details). For A ^> r a reasonable description of the masses could be 
obtained though a consistent picture of D S1 D spectroscopies and decays remained a problem. As the authors noted, 
"the ansatz taken for the potentials V and S may not be as simple as assumed" . The more general form is 
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V SD = XL ■ Si + AtL ■ S 2 + fiS 12 (6) 

only in the particular case of a Coulomb potential need fj, = r|2|. Direct channel couplings (such as to DK and 
D*K thresholds^ ^j) will induce effective potentials that allow the above more general form. Similarly, higher order 
gluon exchange effects in pQCD will also. Indeed, the full spin-dependent structure expected at order a 2 s in QCD 
has been computed 11] and reveals that an additional spin-orbit contribution to the spin-dependent interaction exists 
when quark masses are not equal. When these are incorporated in a constituent quark model there can be significant 
mass shifts leading to a lowered mass for the D s q consistent with the D s o(2317)[l2|- Here we apply this model to the 
recently discovered D s states. 

II. CANONICAL cs STATES 

Predictions of the new model in the D s sector are summarised in Table [I] (the 'high' parameters of Ref . |3| are 
employed) . 



TABLE I: D s Spectrum. 



state 


mass (GeV) 


expt[13] (GeV) 




1.968 


1.968 


D 3 {2 1 So) 
DsftSo) 
D* B (l 3 Si) 


2.637 




3.097 




2.112 


2.112 


D* s (2 3 Si) 
D s *(3 3 Si) 


2.711 


2.688? 


3.153 




A,(l 3 Di) 


2.784 




Dso(l 3 P ) 


2.329 


2.317 


D s0 (2 3 P ) 


2.817 


2.857? 


D a0 (3 3 P ) 


3.219 




D tl (lP) 


2.474 


2.459 


D sl (2P) 


2.940 




D sl (3P) 


3.332 




D' S1 (1P) 


2.526 


2.535 


D' sl (2P) 


2.995 




D' sl (3P) 


3.389 




D s2 (l 3 P 2 ) 


2.577 


2.573 


D s2 {2 3 P 2 ) 


3.041 




D s2 (3 3 P 2 ) 


3.431 





Since the D s j(2690) and D s j(2860) decay to two pseudoscalars, their quantum numbers are J p = + , 1~, 2 + , 
etc. Given the known states 13] and that the energy gap for radial excitation is hundreds of MeV, on almost model 
independent grounds the only possibility for a D s j (2690) is an excited vector. Table |U shows that the D s j(2690) can 
most naturally be identified with the excited vector D*(2S); the D-wave vector is predicted to be somewhat too high 
at 2784 MeV though mixing between these two basis states may be expected. For the D s j (2860), Table [I] indicates 
that this is consistent with the radially excited scalar state D s q(2P). It appears that the D S 2(2P) is too heavy to 
form a viable identification. 



III. DECAY PROPERTIES 

Mass spectra alone are insufficient to classify states. Their production and decay properties also need to be compared 
with model expectations. For example, strong decay widths can be computed with the quark model wavefunctions and 
the strong decay vertex of the 3 Po model. An extensive application of the model to heavy-light mesons is presented 
m Ref. 141. Here we focus on the new BaBar states with the results given in Table ITT1 
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TABLE II: Strong Partial Widths for Candidate D a States. 



state (mass) 


decay mode partial width (MeV) 


D* (2S) (2688) 


DK 


22 




D*K 


78 




D a r] 


1 




D* S T) 


2 




total 


103 


D s0 (2P)(2857) 


DK 


80 




D s -q 


10 




total 


90 


D s2 (2P)(2857) 


DK 


3 




Ds-q 







D*K 


18 




DK* 


12 




total 


33 


D s2 (2P)(3041) 


DK 


1 




D sV 







D*K 


6 




DK* 


47 




D*K* 


76 




total 


130 



A. D sJ (2690) 



The total width of the D* S (2S) agrees very well with the measured width of the D s j(2690) (112 ±37 MeV), lending 
support to this identification. No signal in D s rj is seen or expected, whereas the predicted large D* K partial width 
implies that this state should be visible in this decay mode. The data in D*°(K) — ► D°ir°(K) do not support this 
contention; however, the modes D* + (K) — > D°j(K) and D* + (K) — ► D + ir°(K) show indications of a broad structure 
near 2700 MeVpJ. There is the possibility that l 3 Di mixing with 2 3 Si shift the mass down by 30 MeV to that 
observed and also suppress the D*K mode. For a specific illustration, take the model masses for the 2 3 Si as 2.71GeV 
and l 3 Dx as 2.78 GeV. A simple mixing matrix then yields a solution for the physical states with masses 2.69 GeV 
and its predicted heavy partner at around 2.81 GeV with eigenstates 



|£> s *(2690)) « 
|£>:(2810)) « 

and hence a mixing angle consistent with -0.5 radians 



_(-2|lS) + l|LD)) 



-^=(|LS>+2|LD» 



(7) 




The results of an explicit computation in the 3 Po model are shown in Fig. ^ One sees that a mixing angle of 
approximately -0.5 radians suppresses the D*K decay mode of the low vector (with mass set to 2688 MeV) and 
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produces a total width of approximately 110 MeV, in agreement with the data. The orthogonal state would then 
have a mass around 2.81 GeV and has a significant branching ratio to both DK and D*K, albeit with a broad width, 
greater than 200 MeV. 

In summary, if the D s j(2690) is confirmed as vector resonance, then signals in the D*K channel are expected, 
either in the low lying state (if the mixing is weak) or in a higher vector near 2.8 GeV. 



B. D sJ (2860) 

For the D s j (2860), the D S 2{2P) assignment is further disfavored. At either its model mass of 3041 MeV or at 2860 
MeV the DK mode is radically suppressed, due to the D-wave barrier factor. BaBar see their D s j (2860) signal in 
DK and do not observe it in the D*K decay mode, making the D S 2(2P) assignment unlikely. 

By contrast, the properties of D s j (2860) are consistent with those predicted for the D s o(2P). Within the accuracy 
typical of the 3 Po model for S-wave decays, the total width is in accord with the prediction that the D s q(2P) total 
width is less than that of the excited vectors, and qualitatively in accord with the measured 48 ± 12 MeV. 



C. Radiative Transitions 



The meson assignments made here can be tested further by measuring radiative transitions for these states. Predic- 
tions made with the impulse approximation, with and without nonrelativistic reduction of quark spinors, are presented 
in Table ITTT1 



TABLE III: D a El Radiative Transitions (keV). 



decay mode (mass) 



q 1 (MeV) Non Rel rate Rel Rate 



D s * (25) (2688) -> £> s0 7 
L>*(1D)(2784) -> D s0 7 
D s0 (2P)(2857) £> s *7 
D s2 (2P)(3041) -> £> s * 7 



345 


12.7 


4.6 


428 


116 


82 


648 


13 


0.4 


787 


6.8 


1.9 



IV. PRODUCTION 

The production of the radially excited D s o in B decays can be estimated with ISGW and other formalisms [TH fli^ . 
Since vector and scalar cs states can be produced directly from the W current, the decays B — > D* S {2S)D^ or 
D s o(2P)_D(j) serve as a viable source excited D s states. Computationally, the only differences from ground state D s 
production are kinematics and the excited D s decay constants. 

Production systematics can reveal structural information. For example, the decay B — > Df D~ goes via W emission 
with a rate proportional to Vb c V cs , while W exchange gives rise to B° — > D~K + ~ VbcVud and B° — > DfK~ ~ V c dVb u . 
W exchange is suppressed compared to W emission, thus the expected hierarchy of rates is 



T(B° -> D+D-) > T(B° -► D-K+) > T{B° -► D+K~). (8) 

This suppression of W exchange is confirmed by the data^j| with BR(B° — > DfD~) = (6.5 ± 2.1) x 10~ 3 and 
BR(B° -> DjK+) = (3.1 ± 0.8) x 10~ 5 . The decay to D+K~ has not been observed. 

It is therefore intriguing that the observed rate for B° -> D s (2317) + K- ((4.3 ± 1.5) x 10~ 5 ) is comparable to B° -> 
Dj K + . Assuming accurate data, one must conclude either that this simple reasoning is wrong, the D s (2317)~-ftT + 
mode will be found to be large, or the Z? s (2317) is an unusual state. Searching for the process B° — > D S (2317)~ K + 
is clearly of great interest. 

With the previous warning in mind, we proceed to analyse the production of excited D s states in a variety of 
models. Rates with decay constants set to 1 MeV for D s (2317) and Z? s (2860) production assuming that they are 
simple cs scalar and excited scalar states are presented in Table HVl 

Unfortunately, decay constants cannot be accurately computed at this time. We have evaluated ratios of decay 
constants assuming a simple harmonic oscillator quark model, a Coulomb+linear+hyperfine quark model, and a 
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relativised quark model. The resulting ratio for scalar mesons fall in the range , (2860> ss 0.9 — 1.4. The final 

° ° Jd s (2317) 

estimates of the production of excited scalar D s mesons in B decays are thus 

B -> D S (2860)D 



and 



B -> D S (2317)D 

B -> D S (2860)L>* 
B -> D s (2317)£>* 



= 0.6-1.8 (9) 



= 0.3-0.9. (10) 



Decay Mode 


ISGW 


HQET - Luo & Rosnerfl7] 


Pole[17] 


HQET - ColangeloflS] 


D 3 (2317)L» 
D S (2317).D* 
D S (2860)D 
D s (2860)£>* 


2.78 x 10"'' 
1.06 x 10" 7 
2.09 x 10" 7 
4.57 x 10" 8 


1.95 x 10"' 
8.82 x 10~ 8 
1.72 x 10~ 7 
3.61 x 10" 8 


1.91 x 10"' 
8.79 x 10" 8 
1.66 x 10" 7 
3.55 x 10" 8 


2.24 x 10"' 
1.23 x 10" 7 
1.83 x 10" 7 
4.66 x 10" 8 



TABLE IV: Branching ratios to scalars in different models with decay constants set to 1 MeV 
A similar analysis for vector D* production is presented in Table IV1 



Decay Mode 


ISGW 




HQET - Luo & Rosnerfl7] 


Pole[17] 


HQET - Colangelofl8] 


D*D 


1.97 x 10- 


-7 


1.33 x 10"' 


1.32 x 10"' 


1.57 x 10" ' 


D* S D* 


4.20 x 10" 


7 


3.22 x 10" 7 


3.23 x 10" 7 


4.52 x 10" 7 


£> s (2690)13 


1.01 x 10" 


-7 


8.06 x 10~ 8 


7.77 x 10" 8 


8.79 x 10" 8 


D S (2690)L>* 


4.66 x 10" 


7 


3.55 x 10" 7 


3.49 x 10" 7 


4.65 x 10" 7 



TABLE V: Branching ratios to vectors in different models with decay constants set to 1 MeV 

£ 

Estimating vector decay constant ratios as above yields " p ^. (2 ^ 90> ss 0.7 — 1.1. Finally, predicted ratios of excited 



vector production are 

B -> L> s (2690)£> 
B D* S {2110)D 

and 

B -> D S (2690)D* 
B D*(2110)D* 

We note that Eqn. ^2 agrees well with the earlier prediction of Close and Swanson^ 



0.3-0.7 (11) 



0.5-1.3. (12) 



V. SUMMARY AND CONCLUSIONS 



Given the controversial nature of the D s (2317), establishing a consistent picture of the entire D s spectrum is very 
important. The new states claimed by BaBar can be useful in this regard. We have argued that the six known D s 
and two new states can be described in terms of a constituent quark model with novel spin-dependent interactions. 
Predicted strong decay properties of these states appear to agree with experiment. 

Perhaps the most important tasks at present are (i) discovering the D S 2(2P) state, (ii) searching for resonances in 
D*K and DK* up to 3100 MeV, (iii) analysing the angular dependence of the DK final state in D s j (2860) decay, 
(iv) assessing whether the D s j(2690) appears in the D*K channel, (v) searching for these states in B — * D s jD^*' 
with branching ratios of ~ 10" 3 . 







A. Postscript: Belle discovery 

Subsequent to these calculations, and as this report was being completed, Belief has reported a vector state whose 
mass, width, and possibly production rate and decay characteristics are consistent with our predictions. Specifically, 
their measured mass and total width are M = 2715 ± II+14 MeV and r = H5 ±20132 MeV, in remarkable agreement 
with our predictions. The specific parameters we have used in our analysis are contained within their uncertainties. 
_ Belle [l9| find the new state in B decays, which we have proposed as a likely source. They report Br(B — > 
£>°£>*(2700)) x Br(D* (2700) -> D°K+) = (7.2 ± 1.2±H) ■ 1CT 4 . When compared to the production of the ground 
state vector^ which is Br(B -> D°D* (2112)) = (7.2 ± 2.6) • 10~ 3 , the ratio of production rates in B decay is then 
O(0. l)/Br(D* (2700) — > D°K + ). From our Table II, and assuming flavor symmetry for the strong decay, we predict 
that Br(D* (2690) — > D°K + ) ~ 10% , which within the uncertainties will apply also to the Belle state. Thus the 
absolute production rate, within the large uncertainties, appears to be consistent with that predicted in Section 4. 
If the central value of the Belle mass is a true guide, then a significant branching ratio in D*K would be expected 
(Table II and Fig 1). The orthogonal vector state would then be dominantly ID at 2.78 GeV, but hard to produce 
in B decays. These statements depend on the dynamics underlying 2S-1D mixing, which is poorly understood. It is 
therefore very useful that B decay systematics and the strength of the D*K decay channel in the excited vector D s 
mesons can probe this dynamics. 

Searching for this state in the other advocated modes, and improving the uncertainties, now offers a significant test 
of the dynamics discussed here. 
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